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\ The Jensen theorem is used to derive inequalities for semi-classical tunneling prob- 

abilities for systems involving several degrees of freedom. These Jensen inequalities 
are used to discuss several aspects of sub-barrier heavy-ion fusion reactions. The in- 
equality hinges on general convexity properties of the tunneling coefficient calculated 
with the classical action in the classically forbidden region. 

In complex quantum systems with several active degrees of freedom, one usually 
finds a strong deviation of the tunneling probability from the prediction of simple 
one-dimensional barrier penetration model. Experiments on the fusion of nuclei at 
sub-barrier energies have clearly shown a very large enhancement of the tunneling 
probability when compared to simple one- dimensional barrier model calculations. In 
several important theoretical papers [l|, 0], jl], jl] addressing the tunneling problem 
in systems coupled to a reservoir attempts were made to obtain semi-quantitative, 
albeit important estimates of the effects of the reservoir's degrees of freedom on the 
tunneling dynamics of the subsystem of interest. More detailed numerical calculations 
based on the coupled-channels description of e.g. sub-barrier fusion attempt to give 
q ■ quantitative description within a restricted dimension of the reservoir (the number of 
channels strongly coupled to the entrance channel) 0], 0], (sf, [||. 

Furthermore, very low energy fusion of light nuclei such as 2 H + 2 H has been of 
^ ■ interest over the last two decades in the context of the so-called cold fusion. In this 
endeavour the effect of the environment is important. Recent work on the fusion of 
such light nuclei has indicated that in metals electron screening is enhanced, reducing 
the fusion barrier and accordingly enhancing fusion probability. Of course such light 
ion reactions are of great importance in astrophysics [l0(] and the understanding of the 
effect of the environment on them has been under intensive experimental [llj], (l2| . 



13( and theoretical [14|, [15], [16|, [17[ scrutiny. It would be a useful compliment to 
the above discussion to find general inequalities that involve the tunneling probability 
for a sub-system of the many-degrees-of-freedom system when compared to the sub- 
system alone (with the coupling to the reservoir being averaged). This is the aim of the 
present work. We rely on a general theorem in analysis referred to as the Jensen theorem. 



The Jensen inequality [18J ensures that if F(f(Q) is a functional of a function / (£), 

then (F(f (C)))^ ^ F y(f(C))^j if an d only if F is a convex functional of / within the 

interval in which the average ( )<j is being calculated. One possible way of explicitly 
stating the Jensen inequality is the following: 



£moF{m >F 
£ mo 



if and only if F(f) is a convex functional of / within the interval [a, b], and (j)(C) is any 
positive integrable function. If the convexity turns out to be a concavity, the inequality 
is reversed. 

One immediate consequence of the Jensen inequality is Peierls theorem, which was 
used by Peierls jlof to prove that the canonical partition function, Z((3), defined by 
Z(P) = Tr[exp[—(3H]], is greater or equal than exp [—(3TrH}. 

Using Peierls theorem, R. Johnson and C. Goebel (JG) derived an inequality involving 
the reflection above the barrier in order to assess the effect of breakup on the elastic 



scattering of halo nuclei [19] . That inequality clarified why the reaction cross section 
calculated within the Glauber model is appreciably smaller than that calculated using 
the optical limit of the model, a point first emphasized in [2lj], thus resulting in larger 



radii of halo nuclei. In the following we show that the result of JG [19| and that of 21 
can be considered as a consequence of the Jensen inequality. 

In their above cited work, JG considered the elastic S-matrix element for the I th partial 

wave 

S t (E, C) = exp[2^(£, C)] = exp[/] (2) 
where the phase shift 5i, in the JWKB approximation is given by, 

6t(E,C)= lim | / drXh(r',Q- [ dr'k^ (r'))\ (3) 

™ [J r Jr™ J 

Above, ki(r,Q is the local wave number given by ki(r,() = 
J |§ [E — Vi(r) — F(r)G(£)], kl°\r) is the free particle local wave number, r$ is 

the classical turning point defined by ki(ro, Q = and is the corresponding one for 
the free local wave number. The asymptotic wave number is denoted by k = kf*\r = oo) 
and the mass by pL. 

At high energies, one may expand the local wave number in powers of V ^ +F ^ G ^ 
and retain the leading term. This constitutes the Eikonal approximation considered by 



JG [19j. This approximation gives for the phase shift, 

where the impact parameter b = l+ \/ 2 ■ We consider, as JG, the case where the potential, 
and accordingly the form factor, is purely absorptive (V(r) = — iW(r) and F(r) = )• 
Then the phase shift ^Eikonal becomes pure imaginary and / real. From the Jensen 
inequality and from the fact that 5Eikonai{E,b,Q is a linear function of G(Q we obtain 
the following inequality 



Si(E,£) > exp[-2\S Eikonal {E,b,()\) 



(5) 



which is the result obtained in the work of JG. 

The above inequality only holds for imaginary phase shifts. Clearly the actual 
heavy- ion scattering at intermediate energies involves complex phase shifts, and this 
fact points to an inherent limitation of the work of JG. This limitation is removed if we 
go to very low energies and consider fusion which is dominated by quantum tunneling 
(with real action integral). 

In order to apply the Jensen inequality to fusion reactions, we recall first the expression 
for tunneling probability provided by the coupled-channel treatment in the case of a 
coupling to an oscillator reservoir with zero frequency (sudden approximation), which 
can be cast as a simple average (22| : 

(T,CE7)) f = J dCMOfTiiEMi^+H^irX)) (6) 

where T) [E, Vi (r) + H int (r, £)] is the transmission probability evaluated at energy E with 
an effective potential VJ(r) + H int (r, £), and the wave function, 0o(C)> denotes the ground 
state wave function related to the reservoir coupling. Of course wave functions for excited 
states of the considered reservoir coupling can be used instead. The above equation refers 
to the limit in which the intrinsic energies are small compared to the coupling interaction, 
so that the reservoir Hamiltonian is set equal to zero. 

Using the Kemble [23[ form of the transmission pro bability (25| . which guarantees a 
1/2 transmission at the top of a symmetrical barrier [24|, and takes into account multiple 
reflections inside the barrier to all order if the uniform approximation is used in a path 
integral formulation of tunneling [25], 7} [E, Vi (r) + H int (r, £)] is found to be 

T, [g ,V,( r ) + f ( r ) G (C)] ^ +eTp{a|fi ^. ) + gm[(ri0|} , (7) 
with gi [E, Vi (r) + H int (r, ()) given by 

9l [E, Vi(r) + F (r) G (C)] = J -£ / dry/V^r) + H int (r, Q — E (8) 

V ti J ri (i,o 

where rx(l, C) and r 2 (l, C) are the classical turning points. Here we consider only the case 
where the form factor H int (r, Q is fixed at the position of the maximum of the barrier 
Vi(r), namely H int (r, () = H int (JRj, Q. Although this is a very rough approximation, it 
is a first step in the direction of assessing the effects of the contribution of the coupled 
potential on the transmission coefficient. Bringing the Jensen inequality into the context 
of fusion probability, one can state that 

(9) 

if and only if T [E, V\ (r) + H int (R t , £)] is a convex functional of H int (Ri, ()■ In the 
equation above, (H int (R h £)) f is defined as (H int (R h C)) c = d(\4> (0\ 2 H int (R h Q and 

1 0o (C)| 2 is the square modulus of the normalized ground-state wave function related to 
the reservoir. The same definition holds for the average (Ti(E))^. Hence, it is necessary 
to determine whether the transmission probability is a convex or a concave function 
of H int (RiX) (where H int (Ri,Q is regarded as a simple variable) in order to make a 
comparison of the type of inequality (J9]). The interval [a, b] stands for all possible values 



(T t (E)) c > Ti E, VKr) + (H int (R h ()) t 



that the coordinate related to the oscillator reservoir, (, may assume. 

Let us introduce the quantity w(() = E — H int (Ri,Q , which will be used in our 
calculations in order to make the physical comprehension clearer, that is, w(() will stand 
for the effective energy. Because w(() is a linear function of H int (Ri, () , the sign of the 
second derivative of the tunneling probability 7], with respect to u>(C) determines if 7] is 
a convex functional of the function Hi nt (Ri, () or a concave one: 



d 2 Ti exp [hi (w)] \ r , w . „ 2 . r , (dfi{w) 

(exp [ft, (wjj - 1) (/j (to)) + (exp [Aij (tu)J + r' 



<9w 2 (1 • cxp7/ / (/c) ) i ' - - V dw 

(10) 

in which hl ( W ) = ^frZ] dr^W)^ and fl ( W ) = yf% £™ 

For heavy ions at near-barrier energies, the effective tunneling potential V\ (r) is usually 
approximated by a an inverted parabola, which enables us to follow the Hill- Wheeler 
procedure [26| , in order to obtain a closed form for the Kemble tunneling probability. For 
such heavy ions, the extra degree of freedom, namely the coordinate £, would stand for the 
displacement due to vibrational modes, and the coupled reservoir would be represented 
by an oscillator in this case. Therefore, for such cases, 

Vi (r) = Vnwiir) = Vi(R t ) - ^fjwf(r - Rif (11) 

Hence, 



T2(l,w) 

ri(l,w) 



dr 



Vi(Ri 



Ifjwf (r 



R, 



w 
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dfi H 
dw 



'12) 



This result combined with Eq. (jTOj) yields: 

d 2 T t _ exp [h (w)] 



dw 2 (l + exp[hi{w)]Y 



(exp [hi (w)} - 1) (f t (w)Y > 



(13) 



and finally, 



(T, [w((),V Hwl (r)]) c > T x [(w(C)} c , Vhwi (r) 
or, since (w(()) ( = E - (H int (R h C)) c , 

(Ti [E, V HW i (r) + H mt (R h ()]) c > T t \e, V HWl (r) + (H mt (R h ()) 



(14) 



for all /-partial wave functions and different values of /i. Thus within Hill- Wheeler ap- 
proximation, all systems show enhanced tunneling. 

We now assess the application of the Jensen inequality for the case of low energies. 
By "low" energies, we mean small values of the function w(Q. Here the coupling to the 
reservoir could stand for coupling to the electronic degrees of freedom. The parabolic 
approximation for the potential barrier is not suitable for this case, and we shall use a 
general ion-ion effective interaction, which has the form 



K (r) = V N (r) + — — + ^ 2 (15) 

where Vjy(r) is the nuclear attractive potential. As shown in the Appendix, whatever 
specific analytic form the short-range nuclear interaction, Vjv(r), may take, a potential 
barrier resulting from Vj(r) of Eq.(16) always leads to the following important Jensen 
inequality for very small values of w and/or E: 

(T, VJ (r) + # int C)]) c > Ti [E, V l (r) + (H mt {R h C))J (16) 
where Vj(r) is defined by Eq. ( TT5|) . 

This very general result implies that whatever attractive nuclear potential model 
one may use, the plot of the curve of the transmission probability versus w (£), where 
w (0 = E — H int (Ri, (), is always convex for small values of w, leading to enhanced 
tunneling. 

From the results depicted by Eqs. (15) and by Eq. (17), one is compelled to infer 
that, in general, the tunneling probability 7] would tend to be a convex functional of the 
function w(Q, as seen in Fig.l, in which the tunneling probability was defined by Eq. (J7|), 
and w (() ranges from E m[n to E max , where E min is the minimal energy required for the 
transmission probability to be finite and E max stands for the hight of the tunneling barrier, 
for each considered system. This general property of convexity implies an enhanced 
tunneling or fusion, as experimental data seem to clearly indicate 27j|. 



So far we have concentrated our attention on the I th transmission coefficient. The 
experimental data, on the other hand, are represented by the fusion cross section defined 
by, 

ME) = E( 2/ + w^o = E a ^ ( 17 ) 

" 1=0 1=0 

From Eq. (1T7]) . we see that the dependence of crp(E) on the coupling H int (Ri, £) lies 
only on the terms Ti(E). Then, if it is possible to state that, for instance, Ti(E) is a 
convex functional of H int (Ri,() for all values of the quantum number /, then one can 
also state that <Jf(E) is a convex functional of H int (Ri, £) . The above lends support to 
the general idea that there is an enhancement of the fusion cross section when coupling 
to the degrees of freedom of the reservoir are taken into account, namely 

(<*MO)>c ^ MMC)) C ) (is) 

This is easily seen at deep sub-barrier energies, where in fact the transmission coeffi- 
cient or tunneling probability can be approximated by an exponential, since the action 
in the Kemble formula is small, 

ap{E) = ^E UE) = |S 6XP h9 ° iEl V{r) + Hmt (Rh 0)] (19) 
As shown in the Appendix, the above function is convex in H int (R t , () for u>(C) — > 0, 
and thus its average over ( is greater than that calculated with {H int (Ri, Q)^. Thus, we 
can state, 



(exp [-g (E, V{r) + H int {R h Q)]) > exp [—g (E, V(r) + (H mt (R t , ()),)] (20) 



which represents the very low energy tunneling version of the JG inequality of elastic 
scattering eikonal S-matrix element of halo nuclei. 

In reality, large enhancement in o> has been observed for most heavy- ion fusion 



systems at sub-barrier energies [27|]. Recently, it has been reported that at deep 
sub-barrier energies, this enhancement is reduced 28J (unfortunately, this effect has 
been widely called hindrance, which should not be confused with what we mean by 
hindrance, namely, a cancave behavior of Ti(w(()) as a function of Hi nt (Ri, ()). 



The convexity of the unaveraged tunneling probability for / = can be seen in Fig. 1 
for the systems 64: Ni+ 64 Ni and 16 + 150 Sm. However, plots of the tunneling probability 
as a function of w(Q for very light ions ( 2 H + 2 H, 3 H + 3 H), as shown in Fig. 2, show a 
different behavior. For such light ions, the curve of 7} versus w(() presents an inflection 
point, and thus for high values of w it becomes concave. This result is in contradiction 
with the general analytical result represented by Eq. (fT4l) . where the parabolic potential 
was used to approximate the real potential barrier. That happens possibly because for 
light ions such approximation for the potential barrier is not suitable, as it appears that 
a more accurate approximation that would take into account the highly asymmetrical 
character of the potential curve would be required. A third degree polynomial would be 
a better fit for this purpose, but the analytical treatment becomes extraordinarily more 
complicated. For fusion probabilities involving light ions, the Jensen inequality can only 
be applied within restricted regions of the spectrum of values which w(() may assume. 

In conclusion we have considered some general properties of the tunneling probability 
for systems coupled to a reservoir. Using the Jensen inequality, we have shown that 
within the Kemble/uniform approximation theory of the tunneling probability, the 
average transmission probability is in general larger than that calculated when the 
reservoir degrees of freedom are averaged out at the outset. This has an immediate 
consequence on sub-barrier fusion of heavy ions, where data seem to indicate an 
enhanced tunneling owing to the coupling to the reservoir (coupled channels effects). 
In addition, we have shown that the results obtained by JG [l9| can be generalized by 
using the Jensen inequality. The underlying mathematical dependence of the tunneling 
probability as a function of the reservoir coupling, namely the tunneling probability is 
in general a convex functional of the coupling hamiltonian, permits Jensen inequality 
to be applied to this research field in order to compare two different forms of reaction 
probabilities, both of physical interest. The peculiar behavior presented by the curves 
of transmission probability for light ions indicates that for such systems the effect of 
the coupling to the reservoir on fusion might be different from such effect in heavier 
nuclei, since a change to concavity reverses the Jensen inequality, and the enhancement 
in tunneling becomes a hindrance. 
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APPENDIX 



In this appendix we apply the Jensen inequality for the tunneling probability for very 
small w and/or E, Eq. (17). 

From Eq. (JTU1) . it follows that for small values of w (£), one has 

flTi_ «pfa,(.)] { (fi[w)f + (21) 



dw 2 (1 + exp [hi (w)]y { \ dw 

where fi (w) and hi (w) are defined as in Eq. ({TO]) . From the equation above, we see that 
the sign of will depend exclusively on the term (w)) 2 + y 9 ^q^ j j • We will show 
that such term, considering the potential barrier for fusion reaction with which we are 
dealing (Eq. (|15p ). is always positive when w (() tends to zero. In order to do this, we 

first assume the contrary, namely we suppose that (w)) 2 + ( ^^^ j j ^ 0. Then, 

f d ( 1 M f /2/I /-^(J,™) ^ r 

lim < — — ) > ^ — 1 1 ^ lim {u;/; (w)} =>- 1 ^ lim 




w->0[ dw\fi(w)J) w->0 ' „-* | V H 2 J ri (l,w) y/Vt (r) - w 

Now, let us make 

lim I r (l ' W) dT \ lim \ 1^ dr X lim I dr 

™° \ J ri (i, w) (r) - w J ~ «™ \ 7 ri (,,„) ^ (r) - w J j ^ (r) - w 

in which ri(l,w) < r* < r 2 (l,w). Here r* is chosen to be greater than the distance at 
which the attractive nuclear potential becomes negligible. Hence, for w — > 0, we have 

1 ^ lim { \l%wh \ + lim \ \\%wh \ (22) 

w^o 1 V h 2 I 1 V h 2 f 

where ii = T \ -, dr and J 2 = f r ^ l ' w ^ dr Clearly L is bounded for all values 

Jri(l,w) ^ Vl ^_ w z Jr* ,/v,(r)~w J 

of «; — ► 0, and therefore lim^o <{ \/jji ll 'h \ = 0. That leave us with the inequality 

(23) 

We now turn to the question whether wl 2 is bounded for w — > 0. Performing a change 
of variables, namely y = VTr) — to, one gets for I 2 : 

cfy dV^ 1 (j/ + w) 

Since the point r* is taken to be much greater than the effective nucleus radius, the 
contribution for the total potential V t (r) of the attractive Woods-Saxon potential can 
be neglected within the interval (r* , r 2 (l , w)) . Therefore, in the calculations for I 2 . we 
approximate 




KM « | + ^ (24) 

in which Ci = Z\Zie? and 6*2; = - . Clearly Ci and C21 are non-negative. Here we 
first assume that I 7^ 0, and therefore C21 is strictly positive. From Eq. (12"4"]) . we have 

r _ Cj + VC 1 2 + 4(y + w)C !M 
2 (3/ + w) 

and accordingly 







G 



2/ 



Cx + VC 1 2 + 4(y + «;)C : 



2/ 



(l/ + w) yjy {CI + 4(y + w)G 



21 j 



w 



It is not difficult to prove that J^- > 0. A direct consequence of this fact is that 
limc^o {h (Ci)} ^ h (Ci) , since Ci is positive. Hence 




iy (r*) z 



^ ^2 =>• \/C^ < lim {u>/ 2 (w)} 



2/ 



l«— 



Combining the last result with the inequation ( J23l) we find 

1 ^ 1/?? lim iu>/ 9 Mi ^ 



ft 2 to— V ft 

which implies the absurd result that y/l (I + 1) < 1 since C 2 i = - ■ By assumption 



I 7^ 0, and hence the minimum value for the term a// (/ + 1) is \/2. That leads us to a 

contradiction, and therefore our initial assumption, <Afi {w)) 2 + j ^ can 

not be true. Let us examine now the case where / = 0, which means that the effetive 
potential Vj (r) used in J 2 will be just the Coulomb potential: 



VJr) « — 



With the above potential, the integral J 2 becomes, for w — > 0: 
/ 2 = 



Ci 



to V wr* 



2d 

1 H 3- arctan ^ exp 



«)2 



arccosh 



wr 



ttCi 
3" 



(25) 



Multiplying both sides of Eq. ( 1251) by w and taking the limit w — > 0, we find that 
lim^^o {wh (w)} = 00, and therefore the inequality (l23i) can neither be satisfied for the 
case of the partial wave with I = 0, nor the case / 7^ 0. This proves that the assumption 



we made at the beginning of this section, namely < (fi (w)) + 



dfi(w) 
dw 



^ 0, is false. 



w— 

d 2 T, 



Therefore, recalling Eq. (I2T1) . we have that for small values of w(Q, > 0, which 
implies, for w — > 0, that 



(Ti [E, V { (r) + H mt {R u ())) c >T t E,V t (r) + (H int (R h ()) 



(26) 



where Vi{r) is defined by Eq. ([15 
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FIG. 1: Tunneling probability for I = versus the function w(Q for the systems a) 64 Ni + 
64 Ni, and b) 16 + 150 Sm. The curves show a convex dependence of the tunneling probability 
functional on the function u>(£). 




FIG. 2: Tunneling probability for I = versus the function u>(C)> for the systems a) 2 H+ 2 H, 
and b) 3 H+ 3 H. Both curves show a change in curvature of the tunneling probability functional 
as the function w(£) increases. 



